Vibrations of weakly-coupled nanoparticles 
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The vibrations of a coupled pair of isotropic silver spheres are investigated and compared with the 
vibrations of the single isolated spheres. Situations of both strong coupling and also weak coupling 
are investigated using continuum elasticity and perturbation theory. The numerical calculation of the 
eigenmodes of such dimers is augmented with a symmetry analysis. This checks the convergence and 
applicability of the numerical method and shows how the eigenmodes of the dimer are constructed 
from those of the isolated spheres. The frequencies of the lowest frequency vibrations of such dimers 
are shown to be very sensitive to the strength of the coupling between the spheres. Some of these 
modes can be detected by inelastic light scattering and time-resolved optical measurements which 
provides a convenient way to study the nature of the mechanical coupling in dimers of micro and 
nanoparticles. 
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I. INTRODUCTION 

The ensemble or acoustic vibrations of micro and 
nanoparticles have been studied during the last few 
decades using a variety of spectroscopic and time- 
resolved optical techniques. These vibrations can be used 
for the characterization of the size and shape of the par- 
ticles. Their coupling with the electrons are known to 
contribute to the excitonic dephasingi and even to have 
a universal role in the optical emission of quantum dots? 2 * 
Recently some experimental measurements have focused 
on the acoustic modes of a pair of mechanically coupled 
metallic nanoparticles (NP)^^ or silica microspheres.— 
Measurements also exist for other systems where the me- 
chanical coupling is more complex such as nanocolumns 
made of overlapping spheres^, self-assembled systems for 
which the coherent vibrations of the lattice of spherical 
NPs have been evoked^ and nanopowders where indi- 
vidual NPs are in contact j^t— The goal of the present 
theoretical study of the vibrations of dimers of NPs is to 
improve the understanding of experimental results con- 
cerning dumbbell NPs and also to pave the way to studies 
of a variety of other more complex systems where indi- 
vidual NPs are close enough so that their vibrations can 
be coupled. 

A free homogeneous isotropic sphere has vibrational 
modes which we denote using our notation from a pre- 
vious workii 2 . Briefly, these modes are classified accord- 
ing to four integer quantum numbers. The first is the 
distinction between torsional (T, zero divergence) and 
spheroidal (S, nonzero divergence) modes. Modes are 
next labeled by the usual nonncgative integer angular 
momentum I. There is also an angular momentum z- 
componcnt m. Finally, modes are indexed by integer 
n in order of increasing frequency, starting with 1. We 
thus denote an arbitrary normal mode of a sphere either 

For a silver nanosphere whose diameter is small com- 
pared to the wavelength of light, So and S2 are the 
only Raman active vibrations and So are the only vi- 



brations observed by time-resolved pump-probe experi- 
ments. Since the fundamental vibrations are the main 
features in both kind of experiments, this paper will 
mainly focus on the case of symmetric vibrations orig- 
inating from the fundamental So,o and S2,o- 

II. METHOD 
A. Continuum eigenvibrations and symmetry 

Hathorn et alJ£ used a molecular dynamics model 
to investigate the vibrations of polymer nanoparticle 
dimers. In the present work, we use a continuum elastic 
model which has been shown to be suitable for micro- 
spheres down to rather small nanospheres Moreover, it 
enables a straightforward way to identify the symmetry 
of the eigenvibrations (irreducible representations) and to 
compare with the vibrations of free spheres. We use the 
tools recently presented elsewhere^ 2 - which are suitable 
to calculate the eigenvibrations of arbitrary systems and 
to classify them in terms of symmetry, volume variation 
and projections. The volume integrations required in the 
calculation of the eigenmodes according to the method in- 
troduced by Visschcr et al^ were computed numerically 
by integration only along the z direction by taking ad- 
vantage of the axial symmetry. The wavefunctions were 
expanded in a x l y^z k basis with i + j + k < N. We used 
N = 20 in this work unless stated otherwise. 

The main focus of this work is the case of dimers con- 
sisting of weakly coupled spheres. The eigenmodes of 
any system satisfy two basic rules: they belong to an 
irreducible representation of the point group of interest 
and all the points of the system must oscillate at the 
same frequency for a given eigenmode. These two rules 
let us predict the nature of the eigenmodes in the weak 
coupling regime. 

The point group associated with a dimer NP is Dooh if 
it is made of two identical spheres and Coo if the spheres 
arc different. The irreducible representations arc Ai, A2, 
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Ei, E2, E3, ... for Coo„ and the same with parity (u and 
g) for Dqo/j. The degeneracy with m for systems having 
spherical symmetry is partially lifted for axisymmetric 
systems according to the following rule: S™ modes turn 
into Ai, T™ modes turn into A2 and S^ m ^ and T™ m _^ 
modes turn into E| m | . These rules are obtained by check- 
ing the character table of Coo and remembering the e™^ 
dependence of the displacements. 

The three rigid translations and three rigid rotations 
of the individual spheres can be considered to be vibra- 
tional modes with zero frequency. Conventionally, they 
are not included when enumerating modes. For our pur- 
poses here, it is very helpful to include them. We label 
them as S? m and T° m respectively. This is an exception 
to our convention that modes are indexed with n starting 
from 1. They transform using the rules given before for 
an axisymmetric system. 



B. Perturbation theory 

1. Basic equations 



dynamical equation is 



BiajpUgjp = —Lolm,iU oia 

jp 
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When A is made nonzero, the mode will continuously 
shift to new atomic displacements Ui a and new frequency 
oj. The new dynamical equation is 
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We expand the square of the mode frequency as a 
power series in A as follows: 

uj 2 = lo 2 + Xui 2 + \ 2 u 2 + ... (5) 

and likewise expand the mode displacements in A: 

U la = Uoia + \u Ua + X U2ia + ■■• (6) 

We keep the mode displacements normalized as follows: 
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Consider two elastic objects, where the normal modes 
of vibration are known for each of them individually. In 
this section, we show how a linear perturbation expansion 
can be used to find the frequency shift of modes when 
two such objects are weakly coupled. There are several 
specific situations we have in mind. The first is where 
the two objects are two nanoparticles which are weakly 
coupled together. The second is where a nanoparticle 
is weakly coupled to a substrate. However, this same 
formalism could be applied to a system of three or more 
objects. For clarity, we will always refer to a system of 
two objects below, even though other numbers of objects 
are possible 

The total number of atoms in the system of two objects 
is N. The index i labels the atoms from 1 to N. The 
Cartesian axes are labeled with Greek indices a or j3, 
going from 1 to 3. u.i a is the displacement of atom i 
along the a axis, m, is the mass of atom i. 

Without coupling between the objects, the net force 
on atom i along axis a is 



E B ^P u iP 
iP 



(1) 



where the dynamical matrix of the uncoupled system is 
Biaj/3- When the two objects are coupled, the dynamical 
matrix changes to Ai a jp, where \Ci a jp is the dynami- 
cal matrix due to the coupling alone and A is a scalar 
parameter that we can use for a pcrturbative expansion. 
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Consider a normal mode of the uncoupled system. It 
has frequency ui and atomic displacements u i a . The 



To determine the effect of the coupling term, we sub- 
stitute the expansions in A into Eq. |QJ. Collecting all 
terms linear in A, we obtain 



^2 C iajP u o 3 /3 + ^2 B iajP U ljP 

iP iP 



(8) 



Equation (J5J is now multiplied through by u i a and 
each term is summed over i and a. In addition, substi- 
tution of the expansions in A into Eq. ([7]) and collecting 
linear terms in A tells us that 
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Consequently, we obtain an expression for the mode 
frequency to linear order in A and set A to one. 
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2. Two-point coupling 

We now specialize to the case where the two objects are 
coupled in the simplest possible way. We want to connect 
them by a "spring" . What we mean by "spring" will be 
explained below. In order for the coupling to be weak, 
we restrict its influence to a very small volume fraction 
of both objects as in Ref. H|. We idealize this situation to 
the case of a coupling involving just a single atom on each 
object. Atom a on the first object is coupled through the 
spring to atom b on the second object. Such an arrange- 
ment is only capable of coupling the component of the 
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force which is parallel to the axis of the spring. We will 
suppose that the spring is aligned along the z-axis. Thus, 
atoms a and b share the same x and y coordinates. 

In this situation, there are only four nonzero elements 
of the coupling matrix: C a3a3l C a3b3 , C b3a3 , and C b3b3 . 
Let F a3 denote the z-componcnt of the force on atom a 
from the spring. Then 



and 



F a3 — C a3a3 u a3 + C a3b3 u b3 



F b3 — C b3a3 u a3 + C b3b3 u b3 



(11) 
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The first order perturbation formula for the frequency 
now becomes 



3a3 "L U Q ^C b3b3 + 2u oa3 U ob3 C a 



■Mr.', 



3. Thin cylinder 
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We now consider the case of a "spring" consisting of 
a thin cylinder. If the material has Young's modulus Y 
then the spring constant of this spring is k sp = YA/L 
where A is the cross sectional area of the cylinder and 
L is the length. However, this needs to be modified if 
the frequency of the system is not low compared to the 
internal vibrational modes of the thin cylinder because 
the cylinder does not behave as an ideal massless spring 
in that case. In this section, we analyse the situation 
when the frequency is not necessarily low. 

Atom a is located at the top end of the cylinder, at 
z = L/2. Atom b is located at z = —L/2. u z (z,t) is the 
displacement field inside the cylinder. It has the general 
form 



u z (z,t) = Be* ( **- hrt > + De i <- fc *- wt > 



(14) 



where B and D are constants, k = oj/v and v = \JY~j~p 
is the speed of longitudinal vibrations in a thin rod, 
2872 m/s in silver. 

e zz (z,t) is the strain in the cylinder, given by e zz = 
du z /dz. The stress is a zz (z,t), given by o zz = Ye zz . 
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-Aa zz (L/2,t) 
Aa zz {-L/2,t) 
ut = u z (L/2, t) 
= u z (-L/2,t) 



(15) 
(16) 
(17) 
(18) 



For simplicity, we restrict the remainder of our discus- 
sion to the case where u a3 = —u b3 and F a3 = —F b3 which 
is valid for the symmetric modes we are mainly interested 
in. Note that C a3a 3 - C Q 363 = F a3 /u a3 . Furthermore, 
B + D = and 



u a3 = B(e 



ikL/2 



-ikL /2 \ 



F a3 = —ikYAB(e 



ikL/2 i -ikL/2\ 



C a3 a3 ~ C a3b3 



-ikYA(e lkL / 2 + e' lkL / 2 ) 

„ifcL/2 _ p-ifcL/2 



C a3a3 - C a3b3 = -kYA/tan(kL/2) 
The first order perturbation formula is 
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Finally, we take the continuum limit and apply this to 
the case where the two objects are two identical homoge- 
neous spheres. Each sphere has mass M sp h and volume 
V sp h- In this case, 



pkYA 



M sph tan(fcL/2) 



where we define p as 



(u z (r a )) 2 V sph 
Jsp h u-udV 



(24) 
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Here are some values of p for the north pole of an 
isotropic silver sphere. For the zero frequency translation 
mode along z, p = 1. For the spheroidal mode with m = 
and I equal to 0, 2, 3 and 4, p equals 0.87, 3.26, 6.07 
and 9.18 respectively for the silver spheres considered in 
this work. 



III. RESULTS AND DISCUSSION 

A. Symmetrical dimer 

We consider a dimer made of slightly overlapping iden- 
tical spheres of radius R = 5 nm whose centers are on 
the z axis at z = d/2 and z = —d/2. The perturbation 
approach presented before does not deal with such a cou- 
pling between nanoparticles but its simplicity makes it a 
better starting point to understand how the vibrations 
of a dimer are built. In all this work, we chose to work 
only with an isotropic approximation for silver (mass den- 
sity: 10.5 g/cm 3 , sound speeds: 3747 m/s (longitudinal) 
and 1740 m/s (transverse)). Such an approximation is 
known to be adequate in most cases and in particular in 
the case of multiply twinned particles. Therefore clastic 
anisotropy which has been shown only recently to play a 
significant role for mono-domain gold nanoparticles^ and 
never for silver ones will be ignored here. TableUgives the 
frequencies and irreducible representations of the lowest 
oscillations of a single sphere. Table ILTl gives the assign- 
ments for the lowest frequency modes of a dimer made 
of the same spheres having their center being 9 nm apart 
(d < 2R). 
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v (GHz) 138.5 146.9 202.2 214.0 219.4 281.7 282.2 286.4 319.2 335.1 340.1 347.1 375.8 395.2 396.5 



Tl nl Q 1 ry-il nl cjl rpl 

2 02 Ox ±3 03 O4 I4 



T} Sj s£ Tl 



Ti 



TABLE I. Frequencies of the lowest frequency modes of an isotropic silver sphere having a radius of 5 nm calculated using the 
model by Lamb. 
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decomposition 



rpO 

-Ll.O 

rpO 

L l,±l 
O0 
°1,0 

rpO 

J-l.il 

1 2,0 
T2,±2 
T2,±l 
T2,±2 
S2,±l 
S2,i2 

S 2 ,0 
S2,i2 

S2,il 
T2.il 
B 2,0 

si.o 



eT?, 
e t?, ±1 

© Si.o 
© TS,±i 

e t 2 .o 
e Ti ±2 
© Ti ±1 
© t 2 , ±2 
© s^ ±1 
e s£,± 2 
e si 

© S2,i2 

e s 2 ,±i 

© T 2 ,±i 
© S 2 .o 

© Slo 



So.o © Sj.o 
So.o © Sj.o 
Ss.o © S5.0 



Au (%) 



2.1 
2.2 
1.0 
0.8 
0.0 
0.2 
0.1 
0.1 
0.0 
0.0 
0.0 
0.1 
0.3 
0.2 
0.1 
0.2 

0.1 

0.0 

0.1 

0.0 



TABLE II. Modes of a symmetric dimer made of two over- 
lapping isotropic silver spheres whose radius is 5 nm where 
the distance between the sphere centers is 9 nm. In-phase 
and out-of-phase oscillations of the spheres are labeled with 
© and respectively. The last column shows the decrease of 
the calculated frequencies when increasing from 18 to 20 
so as to illustrate the numerical convergence. 



In the weak coupling regime, the eigenmodes of a dimer 
can be seen as the superposition of one eigenmode for 
each sphere. Weak coupling means that the vibration 
of each sphere is almost unaffected by the presence of 
the other. For a symmetrical dimer, due to the inversion 
symmetry, the eigenmodes of the dimer have to be either 
even or odd which means that the vibrations of both 
spheres have to be identical and either in-phase (identical 
displacements after translating one sphere over the other) 
or out-of-phase (opposite displacements). 

As a first example, we focus on the quadrupolar mode 
o (Ai g for Dqo^). Two eigenmodes of the dimer are 
obtained from this mode: one having the A l£( symme- 
try for which the spheres oscillate in-phase (i — 32 in 





q1 A 

°3.0 1u 



FIG. 1. Odd and even vibrations of the dimer constructed 
from o modes of the free spheres with 1 = 2 (left) and 
£ — 3 (right). The z-axis goes through the centers of both 
spheres and the 3D displacements are obtained by rotation 
around this axis. 



Table [TTJ) and an Ai„ (i = 24) one for out-of-phase os- 
cillations. These modes are illustrated in Fig. [1] Due to 
the coupling, the oscillation is shifted in frequency with 
respect to that of the free spheres. If the vibrations of 
the two spheres result in different displacements for the 
center of the dimer, then the frequency will be increased. 
Conversely, if the two displacements are identical, then 
there will be almost no frequency shift because the dis- 
placement of each sphere is unaffected by the other one. 
By decreasing the coupling between the two spheres (i.e. 
increasing the distance between their centers) the fre- 
quency difference between the two modes of the dimer 
discussed above is reduced as they both evolve towards 
the frequency of the S2 mode of the free spheres. It 
should be noted that close to d = 2R, the various fre- 
quencies do not reach those for a free sphere. This can be 
attributed to the singular nature of this point together 
with the limitation of this continuum elasticity model 
which is unrealistic for very small overlapping volumes. 

Because the parity of the spheroidal vibrations is the 
same as the parity of £, the construction of the odd and 
even vibrations of the dimer is different for odd £ as illus- 
trated in Fig.[TJfor £ = 3 but the resulting even vibration 
of the dimer is always the one having different displace- 
ments for the two spheres at the center of the dimer and 
therefore the one having a larger frequency shift com- 
pared to the frequency of the free spheres. 

The translation along z for a single sphere corresponds 
to a vibration at zero frequency. From this mode, we can 
construct an in-phase oscillation of the dimer which cor- 
responds to the translation along z of the dimer (0 GHz). 
The out-of-phase oscillation corresponds to the lowest 
Aig mode (67.4 GHz). To confirm this assignment we 
checked that the frequency of this Ai a mode tends to 
zero as the coupling between the spheres decreases. In 
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FIG. 2. Frequency of the lowest Ai a mode calculated for two 
silver spheres of radii R — 5 nm connected with a cylinder of 
radius r = R/10 = 0.5 nm and length L (N = 20 full circles, 
N — 18 crosses, N — 16 stars and N = 21 full squares) 
compared to the frequency of two point masses connected 
with a spring having a force constant equal to that of the 
cylinder (line) using Eq. 1241 



order to have a finer control the magnitude of this cou- 
pling, we considered a dumbbell made of spheres con- 
nected by a cylinder whose radius is r = R/10. This 
small dimension ensures that the eigenfrequencies of the 
free cylinder are large which restricts the number of vi- 
brations of the cylinder which can manifest in the low 
frequency range. Figure [5] shows the variation of the fre- 
quency of the lowest frequency K\ g mode as a function 
of the length of the cylinder L. Since the spheres are 
hardly deformed but rather translated during the oscil- 
lation, it is possible to model this oscillation as that of 
two point masses connected by a spring. While a first 
order perturbation theory for a massless spring is possi- 
ble, we present here only a basic derivation valid for this 
particular vibration. The force constant of the massless 
equivalent spring equivalent to the cylinder is k sp = Y-^ 
where Y is the Young's modulus and the cross sectional 
area of the cylinder is A = nr 2 . This expression is valid 
only for L 3> r. The resulting frequency for the dimer is 
then w = y / 2k sp /AI. The very good agreement observed 
in Fig.[2]for L > 4 nm is an indication of both the validity 
of the assignment of this mode and the accuracy of the 
computation of the vibration eigenmodes even for such 
a challenging geometry. This particular oscillation is of 
interest since it was recently experimentally observed^ 
The present work demonstrates that measuring its fre- 
quency is equivalent to determining the force constant 
of the "spring" which connects the two NPs constituting 
the dimer. 

Table [TT] presents a set of low-frequency modes i = 
1 — 12 originating from the translations (S° m ) and rota- 
tions (T° m ) of the spheres. Modes i = 1 — 6 (not shown) 
correspond to the rotations and translations of the dimer 



with zero frequency. These 12 eigenmodes correspond 
to the number of rotations and translations of the two 
spheres. Then modes i = 13 — 32 correspond to combi- 
nations of the S| m or T| m and so on. Again, these 20 
eigenmodes correspond to the number of S2 m and T\ m 
eigenmodes of the two spheres. Modes corresponding to 
the superpositions of the breathing modes (Sj ) are also 
shown. However, since the vibrational density of states 
increases with frequency, even in this relatively low cou- 
pling case this mode mixes significantly with (S5 ). This 
mixing is not apparent in Table [TT] since only the main 
contributions are shown. The same rule regarding the 
shift of the frequencies as discussed for the modes origi- 
nating from before applies for every modes. 

The way the vibrations of a dimer are built is the ana- 
log of the way molecular orbitals are built using the linear 
combination of atomic orbitals (LCAO) method. Such an 
analogy has already been pointed out before in the case of 
acoustic waves in periodic ensemble of spheres in a host 
material^ For example, the Ai u and A± g modes built 
from the breathing mode of the spheres (Sg ) are analo- 
gous to the bonding and anti-bonding molecular orbital 
of H2 built from the Is orbitals. However the reduction 
in the energy of the bonding orbital compared to that of 
the Is atomic orbital does not exist for the Ai„ vibration 
since there is no acoustic equivalence for the screening of 
the charges of the nuclei due to the electrons. 

In order to study the influence of the coupling of the 
vibrations of the two spheres, the variation of the fre- 
quencies of the different vibrations arc plotted in Fig. [3] 
as a function of the distance between the center of the 
spheres. For d = the system consists of a single sphere 
and the frequencies match those reported in Table UJ For 
d very close to 2R, the behavior is singular because for 
d = 2R the two spheres touch at a single point. Moreover, 
for d = 2R~e each sphere feels the other sphere as a mass 
attached to it but it is no longer the case for d = 2R + e. 
So despite the the convergence being good, the interpre- 
tation of the results is not obvious at this point since we 
can expect a non-smooth variation of the frequencies at 
d = 2R. Figure [3] presents the branches corresponding to 
the situation when d < 2R. The lines between the cal- 
culated frequencies connect eigenmodes having the same 
irreducible representation. This was also done in a pre- 
vious work^ The mixings between the branches having 
the same irreducible representations are numerous and 
show that the nature of the vibrations in the case of a 
strong coupling can be quite complex. 



B. Convergence of the numerical method 

The numerical method used in this work has been 
shown to be very accurate for several different geome- 
tries. The question of its reliability for the challenging 
systems in this work is addressed here. To that end, we 
performed the same calculations with N = 18 (Tabic UTJ) 
and N = 16 and N = 21 (Figs. [3] and 2|) to compare with 
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FIG. 3. Variation of the frequencies of a dumbbell made of 
two silver spheres with radii R = 5 nm with varying distance d 
between their centers. The Ai s and Ai„ branches are plotted 
with circles (blue online) and triangles (red online) respec- 
tively. The frequencies of some Sj modes are marked with 
arrows and labeled with t. 



the reference N = 20. This comparison does not rely on 
the mode index i but rather on the irreducible represen- 
tation. For example, ba> for mode i = 10 in Table |TT] 
is the frequency variation for the lowest frequency Ai g 
mode with A" = 20 and N = 18. We varied A" by steps of 
2 because adding or removing one to N does not change 
the frequencies of all the eigenmodes. This is due to the 
inversion symmetry. Changing A" by one adds or removes 
even or odd functions only depending on the parity of N . 
Therefore only the convergence of even or odd modes is 
changed. 

The convergence for the results presented in Table |TT] 
is good. This is due to the use of a rather high value 
for N. The largest variations A.u are obtained for the 
lowest frequency modes. This tendency is confirmed by 
Figs. [2l [4] and [5] except for large values of L for which 
the convergence for the lowest Ai g modes is much bet- 
ter. The worst convergence is for L = 2 nm for which 
the variation between A^ = 20 and AT = 21 in Fig. [5] is 
still quite large (9%). This can be understood when com- 
paring the calculated displacements along the z-axis with 
that expected from the ideal spring model for which both 
spheres are translated in opposite directions and the dis- 
placement varies linearly inside the spring (see Fig. [5| . 
The steep variation of this displacement inside the cylin- 
der together with the flat variation inside the spheres are 
much more difficult to reproduce with power functions 
of z when the length of the spring (L) is small. More- 
over, for L = 2 nm, L/r = 4 so the behavior of the 
cylinder may differ from the spring approximation. The 
variation of the displacement for overlapping spheres is 
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FIG. 4. Frequency of the second lowest Ai 9 mode calculated 
for two silver spheres of radii R = 5 nm connected with a 
narrow cylinder of radius r = _R/10 = 0.5 nm and length L 
(N = 20 full circles, N — 18 crosses, N = 16 stars and N — 21 
full squares) compared to the first order perturbation calcu- 
lation of the coupling of the S^o eigenmodes of the spheres 
through the extensional eigenmode of the thin cylinder (line) 
using Eq. [24] 



not so steep which explains why the convergence is much 
better in that case. Indeed, the spring approximation 
does not hold as can be checked from the displacement 
plotted in Fig. [5] (top). The deformation of each sphere 
is quite significant or in other words the coupling is not 
so weak. Similar comments apply for the other modes 
which are combinations of rotations and translations of 
the spheres with a non-zero frequency except that the 
relevant characteristic value of the cylinder is not related 
to its stretching but rather to torsion and bending. For 
the other modes at higher frequency, the displacements 
inside the spheres are not constant which can be easier to 
reproduce with power functions of x, y and z. Therefore, 
the convergence can be better than for the eigenmodes 
made of the rotations and translations of the two spheres. 

The convergence for the second lowest A\ g vibration 
was also checked using the same method. In that case, 
the calculations are compared with the first order pertur- 
bation results presented before for a thin cylinder. This is 
necessary because the extensional vibration of this cylin- 
der has a frequency similar to the S2 vibration of the 
spheres for L = 2R = 10 nm. The comparison of both 
calculations in Fig. 2] shows the very high accuracy of 
the RUS calculations for long enough cylinders, i.e. when 
the approximations made for the perturbation theory are 
valid. 



C. Non-symmetrical dimer 

In this general case, there are some restrictions on 
which isolated sphere eigenvibrations can be combined 
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FIG. 5. Normalized displacements along the z-axis for the 
lowest frequency Ai s mode of a dimer of overlapping spheres 
with d/R=1.6 (top) and of a cylindrical-necked dumbbell with 
cylinder length L equal to 2 nm (middle) and 10 nm (bot- 
tom). Each sphere has radius 5 nm. The continuous lines 
correspond to the numerical results for N = 20 while the dot- 
ted ones correspond to a constant displacement inside each 
sphere (translation along z) and a linear variation of the dis- 
placement inside the cylinder when it is present. 

to create a vibration cigenmode for the dimer. First, the 
two spheres' vibrations should share the same irreducible 
representation in Coot, (same m and same character (S or 
T) for in = 0) otherwise their superposition would not 
have a well-defined symmetry. Second, their frequencies 
should match because all the points inside the system 
oscillate at the same frequency for an eigcnmode. 

When going from a symmetric dimer to a non- 
symmetric one, the inversion symmetry is lost and the 
two spheres have different sets of vibrational frequencies. 



For example it is possible that the frequency of the Sj 
mode of one sphere might happen to match the frequency 
of the mode of the other. In this case, some modes 
of the dimer as a whole would be hybridizations of these 
two. However, such an exact match of mode frequen- 
cies is very unlikely. Leaving aside such coincidences, in 
the very weak coupling regime the modes of the dimer 
are simply the superposition of one eigenmode for one 
sphere with the other sphere at rest because that is the 
only way to have all the points of the dimer oscillating at 
the same frequency. When the coupling becomes stronger 
the modes can be mixed even more than in the case of 
the symmetric dimer since parity is lost. In that case, 
any modes having the same value of m can mix provided 
their frequencies are close except for m = for which 
torsional and spheroidal modes can not mix. 



IV. CONCLUSION 

We have demonstrated that an exhaustive description 
of the vibrations of a weakly coupled dimer based on 
continuum calculations, perturbation theory and symme- 
try considerations is possible. While the strong coupling 
regime can be investigated using these numerical tools, 
it is of course quite hard to interpret the results since the 
free oscillations of the spheres constituting the dimer are 
strongly modified. The convergence of the calculations 
is discussed in the weak coupling regime thanks to the 
comparison with perturbation calculations and accurate 
results are obtained in most cases. General equations 
were given to extend the perturbation approach to other 
weakly coupled systems. This work paves the way for 
the understanding of the vibrations of more complex sys- 
tems of interest such as chains of overlapping spheres as 
in nanocolumns^ and chains of non-spherical particles. 



* lucien.saviot@u-bourgogne.fr 
t daniel.murray @ubc.ca 

1 T. Takagahara, Phys. Rev. Lett. 71, 3577 (1993). 

2 D. Oron, A. Aharoni, C. de Mello Donega, J. van Rijssel, 
A. Meijerink, and U. Banin, Phys. Rev. Lett. 102, 177402 
(2009). 

3 Y. Li, Q. Zhang, A. V. Nurmikko, and S. Sun, Nano Lett. 
5, 1689 (2005). 

4 A. L. Tchebotareva, M. A. van Dijk, P. V. Ruijgrok, 
V. Fokkema, M. H. S. Hesselberth, M. Lippitz, and M. Or- 
rit, Chem. Phys. Chem. 10, 111 (2009). 

5 T. Dehoux, T. A. Keif, M. Tomoda, O. Matsuda, O. B. 
Wright, K. Ueno, Y. Nishijima, S. Juodkazis, H. Misawa, 
V. Tournat, and V. E. Gusev, Opt. Lett. 34, 3740 (2009). 

6 J. Burgin, P. Langot, A. Arbouet, J. Margueritat, J. Gon- 
zalo, C. N. Afonso, F. Vallee, A. Mlayah, M. D. Rossell, 
and G. Van Tendeloo, Nano Lett. 8, 1296 (2008). 

7 A. Courty, A. Mermet, P. A. Albouy, E. Duval, and M. P. 
Pileni, Nature Materials 4, 395 (2005). 



8 M. Ivanda, M. Buljan, U. V. Desnica, K. Furic, D. Ristic, 
G. C. Righini, and M. Ferrari, Advances in Science and 
Technology 55, 127 (2008). 

9 C. Pighini, D. Aymes, N. Millot, and L. Saviot, J. Nanopar- 
ticle Research 9, 309 (2007). 

10 D. B. Murray, C. H. Netting, L. Saviot, C. Pighini, N. Mil- 
lot, D. Aymes, and H.-L. Liu, J. Nanoelectron. Optoelec- 
tron. 1, 92 (2006). 

11 L. Saviot, C. H. Netting, D. B. Murray, S. Rols, A. Mermet, 
A.-L. Papa, C. Pighini, D. Aymes, and N. Millot, Phys. 
Rev. B 78, 245426 (2008). 

12 L. Saviot and D. B. Murray, Phys. Rev. B 79, 214101 
(2009). 

13 B. C. Hathorn, B. G. Sumpter, D. W. Noid, R. E. Tuzun, 
and C. Yang, J. Phys. Chem. A 106, 9174 (2002). 

14 N. Combe, P.-M. Chassaing, and F. Demangeot, Phys. 
Rev. B 79, 045408 (2009). 

15 W. M. Visscher, A. Migliori, T. M. Bell, and R. A. Reinert, 
J. Acoust. Soc. Am. 90, 2154 (1991). 



8 



H. Portales, N. Goubet, L. Saviot, S. Adichtchev, D. B. 17 M. Kafesaki and E. N. Economou, Phys. Rev. B 52, 13317 
Murray, A. Mermet, E. Duval, and M.-P. Pileni, Proc. (1995). 
Natl. Acad. Sci. U.S.A 105, 14784 (2008). 



